Shared genetic risk between major orofacial cleft phenotypes in an African population.

Nonsyndromic orofacial clefts (NSOFCs) represent a large proportion (70%-80%) of all OFCs. They can be broadly categorized into nonsyndromic cleft lip with or without cleft palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO). Although NSCL/P and NSCPO are considered etiologically distinct, recent evidence suggests the presence of shared genetic risks. Thus, we investigated the genetic overlap between NSCL/P and NSCPO using African genome-wide association study (GWAS) data on NSOFCs. These data consist of 814 NSCL/P, 205 NSCPO cases, and 2159 unrelated controls. We generated common single-nucleotide variants (SNVs) association summary statistics separately for each phenotype (NSCL/P and NSCPO) under an additive genetic model. Subsequently, we employed the pleiotropic analysis under the composite null (PLACO) method to test for genetic overlap. Our analysis identified two loci with genome-wide significance (rs181737795 [p = 2.58E-08] and rs2221169 [p = 4.5E-08]) and one locus with marginal significance (rs187523265 [p = 5.22E-08]). Using mouse transcriptomics data and information from genetic phenotype databases, we identified MDN1, MAP3k7, KMT2A, ARCN1, and VADC2 as top candidate genes for the associated SNVs. These findings enhance our understanding of genetic variants associated with NSOFCs and identify potential candidate genes for further exploration.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway.

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway.

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underly facial shape variation, yet how those in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest even during the late migratory phase results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both Tfap2 members dysregulated numerous midface GRN components involved in midface fusion, patterning, and differentiation. Notably, Alx1/3/4 (Alx) transcript levels are reduced, while ChIP-seq analyses suggest TFAP2 directly and positively regulates Alx gene expression. TFAP2 and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish further implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a mutant zebrafish present abnormal alx3 expression patterns, and the two genes display a genetic interaction in this species. Together, these data demonstrate a critical role for TFAP2 in regulating vertebrate midfacial development in part through ALX transcription factor gene expression.

Transcriptome analyses of murine right and left maxilla-mandibular complex.

OBJECTIVE: The objective of the study was to investigate differential gene expression between murine right and left maxilla-mandibular (MxMn) complexes. SETTING AND SAMPLE POPULATION: Wild-type (WT) C57BL/6 embryonic (E) day 14.5 (n = 3) and 18.5 (n = 3) murine embryos. METHODS: The E14.5 and 18.5 embryos were harvested and hemi-sectioned the MxMn complexes into right and left halves in the mid-sagittal plane. We isolated total RNA using Trizol reagent and further purified using the RNA-easy kit (QIAGEN). We confirmed equal expression of house-keeping genes in right and left halves using RT-PCR and then performed paired-end whole mRNA sequencing in LC Sciences (Houston, TX) followed by differential transcript analyses (>1 or <-1 log fold change; p < .05; q < .05; and FPKM >0.5 in 2/3 samples). The Mouse Genome Informatics and Online Mendelian Inheritance in Man databases as well as gnomAD constraint scores were used to prioritize differentially expressed transcripts. RESULTS: There were 19 upregulated and 19 downregulated transcripts at E14.5 and 8 upregulated and 17 downregulated transcripts at E18.5 time-points. These differentially expressed transcripts were statistically significant and shown to be associated with craniofacial phenotypes in mouse models. These transcripts also have significant gnomAD constraint scores and are enriched in biological processes critical for embryogenesis. CONCLUSIONS: We identified significant differential expression of transcripts between E14.5 and 18.5 murine right and left MxMn complexes. These findings when extrapolated to humans, they may provide a biological basis for facial asymmetry. Further experiments are required to validate these findings in murine models with craniofacial asymmetry.

Computational identification of key transcription factors for embryonic and postnatal Sox2+ dental epithelial stem cell.

While many reptiles can replace their tooth throughout life, human loss the tooth replacement capability after formation of the permanent teeth. It was thought that the difference in tooth regeneration capability depends on the persistence of a specialized dental epithelial structure, the dental lamina that contains dental epithelial stem cells (DESC). Currently, we know very little about DESC such as what genes are expressed or its chromatin accessibility profile. Multiple markers of DESC have been proposed such as Sox2 and Lgr5 . Few single cell RNA-seq experiments have been performed previously, but no obvious DESC cluster was identified in these scRNA-seq datasets, possible due to that the expression level of DESC markers such as Sox2 and Lgr5 is too low or the percentage of DESC is too low in whole tooth. We utilize a mouse line Sox2-GFP to enrich Sox2+ DESC and use Smart-Seq2 protocol and ATAC-seq protocol to generate transcriptome profile and chromatin accessibility profile of P2 Sox2+ DESC. Additionally, we generate transcriptome profile and chromatin accessibility profile of E11.5 Sox2+ dental lamina cells. With transcriptome profile and chromatin accessibility profile, we systematically identify potential key transcription factors for E11.5 Sox2+ cells and P2 Sox2+ cells. We identified transcription factors including Pitx2, Id3, Pitx1, Tbx1, Trp63, Nkx2-3, Grhl3, Dlx2, Runx1, Nfix, Zfp536 , etc potentially formed the core transcriptional regulatory networks of Sox2+ DESC in both embryonic and postnatal stages.

TFAP2 paralogs facilitate chromatin access for MITF at pigmentation and cell proliferation genes.

In developing melanocytes and in melanoma cells, multiple paralogs of the Activating-enhancer-binding Protein 2 family of transcription factors (TFAP2) contribute to expression of genes encoding pigmentation regulators, but their interaction with Microphthalmia transcription factor (MITF), a master regulator of these cells, is unclear. Supporting the model that TFAP2 facilitates MITF's ability to activate expression of pigmentation genes, single-cell seq analysis of zebrafish embryos revealed that pigmentation genes are only expressed in the subset of mitfa-expressing cells that also express tfap2 paralogs. To test this model in SK-MEL-28 melanoma cells we deleted the two TFAP2 paralogs with highest expression, TFAP2A and TFAP2C, creating TFAP2 knockout (TFAP2-KO) cells. We then assessed gene expression, chromatin accessibility, binding of TFAP2A and of MITF, and the chromatin marks H3K27Ac and H3K27Me3 which are characteristic of active enhancers and silenced chromatin, respectively. Integrated analyses of these datasets indicate TFAP2 paralogs directly activate enhancers near genes enriched for roles in pigmentation and proliferation, and directly repress enhancers near genes enriched for roles in cell adhesion. Consistently, compared to WT cells, TFAP2-KO cells proliferate less and adhere to one another more. TFAP2 paralogs and MITF co-operatively activate a subset of enhancers, with the former necessary for MITF binding and chromatin accessibility. By contrast, TFAP2 paralogs and MITF do not appear to co-operatively inhibit enhancers. These studies reveal a mechanism by which TFAP2 profoundly influences the set of genes activated by MITF, and thereby the phenotype of pigment cells and melanoma cells.

AP-2α and AP-2ß cooperatively function in the craniofacial surface ectoderm to regulate chromatin and gene expression dynamics during facial development.

The facial surface ectoderm is essential for normal development of the underlying cranial neural crest cell populations, providing signals that direct appropriate growth, patterning, and morphogenesis. Despite the importance of the ectoderm as a signaling center, the molecular cues and genetic programs implemented within this tissue are understudied. Here, we show that removal of two members of the AP-2 transcription factor family, AP-2α and AP-2ß, within the early embryonic ectoderm of the mouse leads to major alterations in the craniofacial complex. Significantly, there are clefts in both the upper face and mandible, accompanied by fusion of the upper and lower jaws in the hinge region. Comparison of ATAC-seq and RNA-seq analyses between controls and mutants revealed significant changes in chromatin accessibility and gene expression centered on multiple AP-2 binding motifs associated with enhancer elements within these ectodermal lineages. In particular, loss of these AP-2 proteins affects both skin differentiation as well as multiple signaling pathways, most notably the WNT pathway. We also determined that the mutant clefting phenotypes that correlated with reduced WNT signaling could be rescued by Wnt1 ligand overexpression in the ectoderm. Collectively, these findings highlight a conserved ancestral function for AP-2 transcription factors in ectodermal development and signaling, and provide a framework from which to understand the gene regulatory network operating within this tissue that directs vertebrate craniofacial development.

The Skull's Girder: A Brief Review of the Cranial Base.

The cranial base is a multifunctional bony platform within the core of the cranium, spanning rostral to caudal ends. This structure provides support for the brain and skull vault above, serves as a link between the head and the vertebral column below, and seamlessly integrates with the facial skeleton at its rostral end. Unique from the majority of the cranial skeleton, the cranial base develops from a cartilage intermediate-the chondrocranium-through the process of endochondral ossification. Owing to the intimate association of the cranial base with nearly all aspects of the head, congenital birth defects impacting these structures often coincide with anomalies of the cranial base. Despite this critical importance, studies investigating the genetic control of cranial base development and associated disorders lags in comparison to other craniofacial structures. Here, we highlight and review developmental and genetic aspects of the cranial base, including its transition from cartilage to bone, dual embryological origins, and vignettes of transcription factors controlling its formation.

Anomalous incisor morphology indicates tissue-specific roles for Tfap2a and Tfap2b in tooth development.

Mice possess two types of teeth that differ in their cusp patterns; incisors have one cusp and molars have multiple cusps. The patterning of these two types of teeth relies on fine-tuning of the reciprocal molecular signaling between dental epithelial and mesenchymal tissues during embryonic development. The AP-2 transcription factors, particularly Tfap2a and Tfap2b, are essential components of such epithelial-mesenchymal signaling interactions that coordinate craniofacial development in mice and other vertebrates, but little is known about their roles in the regulation of tooth development and shape. Here we demonstrate that incisors and molars differ in their temporal and spatial expression of Tfap2a and Tfap2b. At the bud stage, Tfap2a is expressed in both the epithelium and mesenchyme of the incisors and molars, but Tfap2b expression is restricted to the molar mesenchyme, only later appearing in the incisor epithelium. Tissue-specific deletions show that loss of the epithelial domain of Tfap2a and Tfap2b affects the number and spatial arrangement of the incisors, notably resulting in duplicated lower incisors. In contrast, deletion of these two genes in the mesenchymal domain has little effect on tooth development. Collectively these results implicate epithelial expression of Tfap2a and Tfap2b in regulating the extent of the dental lamina associated with patterning the incisors and suggest that these genes contribute to morphological differences between anterior (incisor) and posterior (molar) teeth within the mammalian dentition.

Finding MEMO-Emerging Evidence for MEMO1's Function in Development and Disease.

Although conserved throughout animal kingdoms, the protein encoded by the gene Mediator of ERBB2 Driven Cell Motility 1 or MEMO1, has only recently come into focus. True to its namesake, MEMO1 first emerged from a proteomic screen of molecules bound to the ERBB2 receptor and was found to be necessary for efficient cell migration upon receptor activation. While initially placed within the context of breast cancer metastasis-a pathological state that has provided tremendous insight into MEMO1's cellular roles-MEMO1's function has since expanded to encompass additional cancer cell types, developmental processes during embryogenesis and homeostatic regulation of adult organ systems. Owing to MEMO1's deep conservation, a variety of model organisms have been amenable to uncovering biological facets of this multipurpose protein; facets ranging from the cellular (e.g., receptor signaling, cytoskeletal regulation, redox flux) to the organismal (e.g., mineralization and mineral homeostasis, neuro/gliogenesis, vasculogenesis) level. Although these facets emerge at the intersection of numerous biological and human disease processes, how and if they are interconnected remains to be resolved. Here, we review our current understanding of this 'enigmatic' molecule, its role in development and disease and open questions emerging from these previous studies.

Differential 3' Processing of Specific Transcripts Expands Regulatory and Protein Diversity Across Neuronal Cell Types.

Alternative polyadenylation (APA) regulates mRNA translation, stability, and protein localization. However, it is unclear to what extent APA regulates these processes uniquely in specific cell types. Using a new technique, cTag-PAPERCLIP, we discovered significant differences in APA between the principal types of mouse cerebellar neurons, the Purkinje and granule cells, as well as between proliferating and differentiated granule cells. Transcripts that differed in APA in these comparisons were enriched in key neuronal functions and many differed in coding sequence in addition to 3'UTR length. We characterize Memo1, a transcript that shifted from expressing a short 3'UTR isoform to a longer one during granule cell differentiation. We show that Memo1 regulates granule cell precursor proliferation and that its long 3'UTR isoform is targeted by miR-124, contributing to its downregulation during development. Our findings provide insight into roles for APA in specific cell types and establish a platform for further functional studies.

AP-2α and AP-2ß cooperatively orchestrate homeobox gene expression during branchial arch patterning.

The evolution of a hinged moveable jaw with variable morphology is considered a major factor behind the successful expansion of the vertebrates. DLX homeobox transcription factors are crucial for establishing the positional code that patterns the mandible, maxilla and intervening hinge domain, but how the genes encoding these proteins are regulated remains unclear. Herein, we demonstrate that the concerted action of the AP-2α and AP-2ß transcription factors within the mouse neural crest is essential for jaw patterning. In the absence of these two proteins, the hinge domain is lost and there are alterations in the size and patterning of the jaws correlating with dysregulation of homeobox gene expression, with reduced levels of Emx, Msx and Dlx paralogs accompanied by an expansion of Six1 expression. Moreover, detailed analysis of morphological features and gene expression changes indicate significant overlap with various compound Dlx gene mutants. Together, these findings reveal that the AP-2 genes have a major function in mammalian neural crest development, influencing patterning of the craniofacial skeleton via the DLX code, an effect that has implications for vertebrate facial evolution, as well as for human craniofacial disorders.

Beyond MITF: Multiple transcription factors directly regulate the cellular phenotype in melanocytes and melanoma.

MITF governs multiple steps in the development of melanocytes, including specification from neural crest, growth, survival, and terminal differentiation. In addition, the level of MITF activity determines the phenotype adopted by melanoma cells, whether invasive, proliferative, or differentiated. However, MITF does not act alone. Here, we review literature on the transcription factors that co-regulate MITF-dependent genes. ChIP-seq studies have indicated that the transcription factors SOX10, YY1, and TFAP2A co-occupy subsets of regulatory elements bound by MITF in melanocytes. Analyses at single loci also support roles for LEF1, RB1, IRF4, and PAX3 acting in combination with MITF, while sequence motif analyses suggest that additional transcription factors colocalize with MITF at many melanocyte-specific regulatory elements. However, the precise biochemical functions of each of these MITF collaborators and their contributions to gene expression remain to be elucidated. Analogous to the transcriptional networks in morphogen-patterned tissues during embryogenesis, we anticipate that the level of MITF activity is controlled not only by the concentration of activated MITF, but also by additional transcription factors that either quantitatively or qualitatively influence the expression of MITF-target genes.

TFAP2 paralogs regulate melanocyte differentiation in parallel with MITF.

Mutations in the gene encoding transcription factor TFAP2A result in pigmentation anomalies in model organisms and premature hair graying in humans. However, the pleiotropic functions of TFAP2A and its redundantly-acting paralogs have made the precise contribution of TFAP2-type activity to melanocyte differentiation unclear. Defining this contribution may help to explain why TFAP2A expression is reduced in advanced-stage melanoma compared to benign nevi. To identify genes with TFAP2A-dependent expression in melanocytes, we profile zebrafish tissue and mouse melanocytes deficient in Tfap2a, and find that expression of a small subset of genes underlying pigmentation phenotypes is TFAP2A-dependent, including Dct, Mc1r, Mlph, and Pmel. We then conduct TFAP2A ChIP-seq in mouse and human melanocytes and find that a much larger subset of pigmentation genes is associated with active regulatory elements bound by TFAP2A. These elements are also frequently bound by MITF, which is considered the "master regulator" of melanocyte development. For example, the promoter of TRPM1 is bound by both TFAP2A and MITF, and we show that the activity of a minimal TRPM1 promoter is lost upon deletion of the TFAP2A binding sites. However, the expression of Trpm1 is not TFAP2A-dependent, implying that additional TFAP2 paralogs function redundantly to drive melanocyte differentiation, which is consistent with previous results from zebrafish. Paralogs Tfap2a and Tfap2b are both expressed in mouse melanocytes, and we show that mouse embryos with Wnt1-Cre-mediated deletion of Tfap2a and Tfap2b in the neural crest almost completely lack melanocytes but retain neural crest-derived sensory ganglia. These results suggest that TFAP2 paralogs, like MITF, are also necessary for induction of the melanocyte lineage. Finally, we observe a genetic interaction between tfap2a and mitfa in zebrafish, but find that artificially elevating expression of tfap2a does not increase levels of melanin in mitfa hypomorphic or loss-of-function mutants. Collectively, these results show that TFAP2 paralogs, operating alongside lineage-specific transcription factors such as MITF, directly regulate effectors of terminal differentiation in melanocytes. In addition, they suggest that TFAP2A activity, like MITF activity, has the potential to modulate the phenotype of melanoma cells.

MEMO1 drives cranial endochondral ossification and palatogenesis.

The cranial base is a component of the neurocranium and has a central role in the structural integration of the face, brain and vertebral column. Consequently, alteration in the shape of the human cranial base has been intimately linked with primate evolution and defective development is associated with numerous human facial abnormalities. Here we describe a novel recessive mutant mouse strain that presented with a domed head and fully penetrant cleft secondary palate coupled with defects in the formation of the underlying cranial base. Mapping and non-complementation studies revealed a specific mutation in Memo1 - a gene originally associated with cell migration. Expression analysis of Memo1 identified robust expression in the perichondrium and periosteum of the developing cranial base, but only modest expression in the palatal shelves. Fittingly, although the palatal shelves failed to elevate in Memo1 mutants, expression changes were modest within the shelves themselves. In contrast, the cranial base, which forms via endochondral ossification had major reductions in the expression of genes responsible for bone formation, notably matrix metalloproteinases and markers of the osteoblast lineage, mirrored by an increase in markers of cartilage and extracellular matrix development. Concomitant with these changes, mutant cranial bases showed an increased zone of hypertrophic chondrocytes accompanied by a reduction in both vascular invasion and mineralization. Finally, neural crest cell-specific deletion of Memo1 caused a failure of anterior cranial base ossification indicating a cell autonomous role for MEMO1 in the development of these neural crest cell derived structures. However, palate formation was largely normal in these conditional mutants, suggesting a non-autonomous role for MEMO1 in palatal closure. Overall, these findings assign a new function to MEMO1 in driving endochondral ossification in the cranium, and also link abnormal development of the cranial base with more widespread effects on craniofacial shape relevant to human craniofacial dysmorphology.

The old and new face of craniofacial research: How animal models inform human craniofacial genetic and clinical data.

The craniofacial skeletal structures that comprise the human head develop from multiple tissues that converge to form the bones and cartilage of the face. Because of their complex development and morphogenesis, many human birth defects arise due to disruptions in these cellular populations. Thus, determining how these structures normally develop is vital if we are to gain a deeper understanding of craniofacial birth defects and devise treatment and prevention options. In this review, we will focus on how animal model systems have been used historically and in an ongoing context to enhance our understanding of human craniofacial development. We do this by first highlighting "animal to man" approaches; that is, how animal models are being utilized to understand fundamental mechanisms of craniofacial development. We discuss emerging technologies, including high throughput sequencing and genome editing, and new animal repository resources, and how their application can revolutionize the future of animal models in craniofacial research. Secondly, we highlight "man to animal" approaches, including the current use of animal models to test the function of candidate human disease variants. Specifically, we outline a common workflow deployed after discovery of a potentially disease causing variant based on a select set of recent examples in which human mutations are investigated in vivo using animal models. Collectively, these topics will provide a pipeline for the use of animal models in understanding human craniofacial development and disease for clinical geneticist and basic researchers alike.

New Functional Signatures for Understanding Melanoma Biology from Tumor Cell Lineage-Specific Analysis.

Molecular signatures specific to particular tumor types are required to design treatments for resistant tumors. However, it remains unclear whether tumors and corresponding cell lines used for drug development share such signatures. We developed similarity core analysis (SCA), a universal and unsupervised computational framework for extracting core molecular features common to tumors and cell lines. We applied SCA to mRNA/miRNA expression data from various sources, comparing melanoma cell lines and metastases. The signature obtained was associated with phenotypic characteristics in vitro, and the core genes CAPN3 and TRIM63 were implicated in melanoma cell migration/invasion. About 90% of the melanoma signature genes belong to an intrinsic network of transcription factors governing neural development (TFAP2A, DLX2, ALX1, MITF, PAX3, SOX10, LEF1, and GAS7) and miRNAs (211-5p, 221-3p, and 10a-5p). The SCA signature effectively discriminated between two subpopulations of melanoma patients differing in overall survival, and classified MEKi/BRAFi-resistant and -sensitive melanoma cell lines.

Transcription factor MITF and remodeller BRG1 define chromatin organisation at regulatory elements in melanoma cells.

Microphthalmia-associated transcription factor (MITF) is the master regulator of the melanocyte lineage. To understand how MITF regulates transcription, we used tandem affinity purification and mass spectrometry to define a comprehensive MITF interactome identifying novel cofactors involved in transcription, DNA replication and repair, and chromatin organisation. We show that MITF interacts with a PBAF chromatin remodelling complex comprising BRG1 and CHD7. BRG1 is essential for melanoma cell proliferation in vitro and for normal melanocyte development in vivo. MITF and SOX10 actively recruit BRG1 to a set of MITF-associated regulatory elements (MAREs) at active enhancers. Combinations of MITF, SOX10, TFAP2A, and YY1 bind between two BRG1-occupied nucleosomes thus defining both a signature of transcription factors essential for the melanocyte lineage and a specific chromatin organisation of the regulatory elements they occupy. BRG1 also regulates the dynamics of MITF genomic occupancy. MITF-BRG1 interplay thus plays an essential role in transcription regulation in melanoma.

A polymorphism in IRF4 affects human pigmentation through a tyrosinase-dependent MITF/TFAP2A pathway.

Sequence polymorphisms linked to human diseases and phenotypes in genome-wide association studies often affect noncoding regions. A SNP within an intron of the gene encoding Interferon Regulatory Factor 4 (IRF4), a transcription factor with no known role in melanocyte biology, is strongly associated with sensitivity of skin to sun exposure, freckles, blue eyes, and brown hair color. Here, we demonstrate that this SNP lies within an enhancer of IRF4 transcription in melanocytes. The allele associated with this pigmentation phenotype impairs binding of the TFAP2A transcription factor that, together with the melanocyte master regulator MITF, regulates activity of the enhancer. Assays in zebrafish and mice reveal that IRF4 cooperates with MITF to activate expression of Tyrosinase (TYR), an essential enzyme in melanin synthesis. Our findings provide a clear example of a noncoding polymorphism that affects a phenotype by modulating a developmental gene regulatory network.

Gene regulatory evolution and the origin of macroevolutionary novelties: insights from the neural crest.

The appearance of novel anatomic structures during evolution is driven by changes to the networks of transcription factors, signaling pathways, and downstream effector genes controlling development. The nature of the changes to these developmental gene regulatory networks (GRNs) is poorly understood. A striking test case is the evolution of the GRN controlling development of the neural crest (NC). NC cells emerge from the neural plate border (NPB) and contribute to multiple adult structures. While all chordates have a NPB, only in vertebrates do NPB cells express all the genes constituting the neural crest GRN (NC-GRN). Interestingly, invertebrate chordates express orthologs of NC-GRN components in other tissues, revealing that during vertebrate evolution new regulatory connections emerged between transcription factors primitively expressed in the NPB and genes primitively expressed in other tissues. Such interactions could have evolved by two mechanisms. First, transcription factors primitively expressed in the NPB may have evolved new DNA and/or cofactor binding properties (protein neofunctionalization). Alternately, cis-regulatory elements driving NPB expression may have evolved near genes primitively expressed in other tissues (cis-regulatory neofunctionalization). Here we discuss how gene duplication can, in principle, promote either form of neofunctionalization. We review recent published examples of interspecies gene-swap, or regulatory-element-swap, experiments that test both models. Such experiments have yielded little evidence to support the importance of protein neofunctionalization in the emergence of the NC-GRN, but do support the importance of novel cis-regulatory elements in this process. The NC-GRN is an excellent model for the study of gene regulatory and macroevolutionary innovation.